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ABSTRACT 

We show that gas in a large fraction of low mass dark matter halos may form 
Toomrc stable disks, if angular momentum is conserved when the gas contracts. Such 
halos would be stable to star formation and therefore remain dark. This may poten- 
tially explain the discrepancy between the predicted and observed number of dwarf 
satellites in the Local Group, as well as the deviation between the predicted and the 
observed faint end slope of the luminosity function. The above mechanism does not 
require a strong variation of the baryon fraction with the virial mass of the dark halo. 
We show that model fits to rotation curves are also consistent with this hypothesis: 
none of the observed galaxies lie in the region of parameter space forbidden by the 
Toomre stability criterion. 
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1 INTRODUCTION 

There is now growing evidence that the observed number of 
dwarf galaxies surrounding the Milky Way or M31 does not 
agree with the predictions of high-resolution N - body simu — 



ago (e.g., White & Rees (1978)). This process expels the gas 
from the halo and therefore suppresses star formation. The 
accretion and cooling of gas can also be suppressed in the 
presence of a strong photo-ionizing background. This route 



has also been known for a number of years now (Doroshke 



lations Hhl trie currently favorer! Hold Dark Matter cosmo- 



vich, Zel'Dovich & Novikov 1967 ) and has been investigated 



l ogical model (flat, low matter-density Universe - LCDM) 



(Klypiq-et- ol. lOOO i-M eoro ot al. 1999). Given the 



of the durrent LCDM paradi gm in reproduci n g many of the 



large- s cale observables (e . g 



recently in the context of semi-analytic models of galaxy 
form a tion (|Bullock ct al. 200 j ; |Chiu et al. 200l[ |Somerville 
200 1| ; [Benson et al. 2001a|) - 



Jaffe et al (2001 



(2001); Verde et al. (2001); Lahav et al. (200l|)), it seems nat 



Peacock et a] 



ural to try to make changes to the existing paradigm without 
modifying its large-scale properties. There are two routes in 
this direction: one is to prevent the formation of small dark 
matter halos, the second is to hide them by suppressing star 
forma.tipn. Some modifications belonging to the first group 



A closely related problem is the fact that the observed 
f aint end of the lum inosity function has a shallower slope 
(Blanton et al. 2001) than predicted from high-resolution 
LCDM simulations, if one simply assume s a constant mass- 
to-light ratio (e.g., Jenkins et al. (2001)). Stellar feedback 
is again often used in semi-analytic models to suppress star 
formation in halos with shallow potential wells (e.g 



consist lof changing the nature of the dark matter itself to 



have different phys ical prop erties: self-interacting ( gpergel 
& Steinhardt 2000 ), warm (Colomb i, Dodelson & Widrow 



1996 Bode, Ostriker fc Turok 2001|) or other variant s ( e.g 



et al. (2000)) 



Cole 



Goodman (200C); Hu, Barkana & Gruzinov (2000) 



Ccn 



(2001)). It has also been proposed that the solution to the 



problem may lie within the nature of the initial conditions 
and that some modification to the shape of the inflaton po- 



tential may suppress small scale structure ( Kamionkowski 



Liddle 2000 ), although this solution faces some problems in 
reproducing the observed amount of structure in the Lyman- 
q forest. On the other hand, the solution of suppressing star 
formation via supernova feedback was suggested a long time 



All of the above models seek to systematically increase 
the mass-to-light ratio in small halos, by depleting the cold 
gas fraction available for star formation. In this paper we 
show that such drastic expulsion of gas is, perhaps, un- 
necessary: provided that the gas conserves angular momen- 
tum during collapse, a large fracti o n of low mass halos 
will be Toomre stable (Toomre 1964; Kcnnicutt 1989| ), i.e., 
Q — E c /E > 1 (where E is the disk surface mass density 
and E c is the critical surface density to trigger gravitational 
instability and thus start formation) . If the baryon fraction 
(fd) in disks is fa ~ 0.5f2;,/f2 m , all disks in halos with masses 
M < 10 9 Mq will be Toomre stable and fail to form stars. 
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The "missing" dwarf galaxies predicted to be around the 
Milky Way and M31 may fall in this category. We also show 
that, if a small fraction of baryons initially present in the 
dark matter halo settle in a disk, as some observational evi- 



den ce suggests ( Jimenez, Verde fc Oh 2002[ Guzik fe Scljak 



2001|; Mm den Bosch, Burkert fc Swatcrs 2001^ |Burkcrt 200(f 
Fukugita, Hogan & Peebles 1998), even more massive halos 



have some probability of remaining dark, provided they have 
high spin. This could potentially reconcile the observed lu- 
minosity function with the mass function predicted in high- 
resolution LCDM N-body simulations. We show that the 
predictions of this model are consistent wi th our previous 
analysi s of t he rotation curves of galaxies (Jimenez, Verde 
& Oh ( |2002| ), hereafter JVO02): none of the galaxies lie in 
the region of parameter space forbidden by the Toomre in- 
stability criterion. 

Throughout this paper we assume a LCDM cosmology 
given by:(n m , fi A , fife, h, agh' 1 ) = (0.3, 0.7, 0.039, 0.7, 1.0). 



2 THE MASS DISTRIBUTION OF DARK 
GALAXIES 

We model spiral galaxies as expon ential disks embedded i n 
a cold dark matter h al o. Following Dalcanton et al. (1997 ) ; 
Jimenez ct al. (1997|); |Mo et al. (1998|) we assume a NFW 
(Navarro, Frenk & White 1997) profile p cx [(r/r c )(l + 
r/r c ) ]~ where r c is the break radius and the surface mass 
density of the disk is given by E = Eo exp[— r/Rd], where 



Rd is the scale length of the disk. 

For a baryonic disk to be locally gravitationally unsta- 
ble, despite the stabilizing effects of tidal shears and pressure 
forces, we require the Toomre parameter Q < 1, where: 



Q = 



(i) 



ttGE' 

E is the disk surface mass density, k = 1.41(V/r)(l + 
dlny/dlnr) 1 / 2 is the epicyclic frequency and c 3 is the gas 
sound speed. 

Regions in local disk galaxies where Q > 1 are observa- 
tionally associated with very little star formation, indicating 
that the Toomre criterion is obeyed remar kably well, if a ga s 
sound speed of c s ~ 6kms -1 is assumed (Kennicutt 1989). 
Ferguson et al. (1996) noted that for one galaxy (NGC 6946) 
the agreement with the Toomre criterion was not perfect: 
there was no truncation of star formation for Q > 1.6 if the 
observed vertical velocity dispersion was adopted. However, 
agreement with the Toomre criterion could be obtained if 
the v elocity dispersion wa s assumed to be constant with ra- 
dius. Wong & Blitz (2001) have measured the value of Q in 
the inner part of rotation curves for a sample of 8 galaxies. 
For half of the sample the predicted value from Toomre the- 
ory for the abundance of gas is in excellent agreement with 
observations. In the other cases the agreement is not as re- 
markable, hut still there is no evidence for large deviations 



from tl c Toomre stability predictions. 



Jimenez et al. (1997) developed a disk model to study 
the dependence of Q on mass, radius and spin parameter 
of the dark matter halo and found that for a certain range 
of the above parameters, some disks would be dark. An im- 
portant ingredient in their study was the assumption of a 
constant baryon fract ion, namely they ad opted the nucle- 
osynthesis value from Walker et al. (1991 ). Using this con- 
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Figure 1. The present day velocity function. Lines labeled by 
Sloan and SRSS2 are obtained from the SDSS and the SRSS2 
surveys respectively. The solid line label ed = 0.06 is the theo 



retical velocity function derived from the Sheth & Tormen (1999) 
mass function using a constant disk baryon fraction = 0.06. 
The line labeled NDG (no dark galaxies) is obtained assuming our 
ansatsz for the dependence of f d on mass (see figure 2). Finally 
the thick solid line labeled by DG (dark galaxies) is obtained con- 
sidering that only disks that are Toomre unstable will be able to 
form stars and thus be visible (see text for more details). 



stant disk fraction for all galaxies, they concluded that for 
an Einstein-de Sitter Universe the minimum dark matter 
halo mass for a visible disk galaxy is about 10 9 Mq . On the 
other hand, for a low matter density universe (e.g., LCDM), 
this minimum mass drops substantially, by about two or- 
ders of magnitude, and thus the fraction of dark halos also 
falls drastically. The reason for this is that for a low density 
Universe with baryon density fixed by nucleosynthesis, the 
baryons constitute a higher fraction of the halo mass. 

It seems unlikely that all baryons will be able to cool 
and form stars. Recent studies from galaxy-galaxy lensing 
from the Sloan Digital Sky Survey (SDSS) show that for 
L* galaxies the baryon fraction is abo ut half of the nucle- 
osynthesis value (Guzik & Seljak 2001). On the other hand, 
theoretical considerations suggest that the baryon fraction 
for galaxies, either less massive or more massive than L», 
is likely to be much sma ller due to feedback a nd inefficient 
cooling respectively (e.g. ra,n den Bosch (2001); 



Benson et al 



(2001b)). This implies that a large fraction of disks should 
be Toomre stable. Below we quantify how this affects the 
predicted number density of observable galaxies. 



2.1 Luminosity function 

CDM models typically predict many more galaxies at the 
faint end of the luminosity function than are actually ob- 



The abundance of dark galaxies 3 



served. This can most easily be seen by comparing the 
galaxy circular velocity function with theoretical predic- 
tions. The galaxy velocity function is obtained by combin- 
ing the observed luminosity function with empirically deter- 
mined luminosity-velocity relations (Tully- Fisher relation 



for spirals, and Faber- Jackson for ellipticals). Gonzalez et al 
(200(|~Jind that treating the entire galaxy population as spi- 
rals does not significantly alter the derived velocity function. 
In fig. [l], we plot a representative velocity function fr om their 
paper, derived from the SSRS2 (|Marzke et al. 1998b B band 
survey and the Yasuda et al. (1997 ) Tully-Fisher relation. 
We also construct a new velocity function from the Sloan 
Digital Sky Survey (SPS S), which has publis hed luminos- 
ity functions in 5 bands ( Blanton et al. 2001 ). We choose 
to use the R band, as the R band Tully-Fisher relation is 
tighter and has less s catter than the Tully -Fisher relation at 
shorter wavelengths. Blanton et al. (2001) perfor m the color 
transfo rma tions from SDSS magnitudes to LCRS (Bhectman 
et al. 1996) R-band magnitudes (see their Table 3); we per- 
form the color transformations and corrections for internal 
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extinction as in |Gonzalcz et al. (200C ) to the -Rcour tcau mag- 
nitudes used by pourteau (1997 ). We then use the Courteau 
(1997^~^ C ourtcau band Tully-Fisher relation to obtain the ve- 
locity function. Note that the Sloan velocity function has a 
significantly higher normalization than the SSRS2 velocity 
function. This is due to the use of Petrosian magnitudes 




for SDSb, wnicn allows detection ot galaxies ot sigmncantly 
lower sAfface brightness. The luminosity density of galaxies 
in the SDSS is therefore significantly higher than other red- 
shift surveys. The SDSS luminosity function is very similar 
to the luminosity functions derived by other surveys such as 
the Las Campanas Redshift Survey (LCRS) and Two Degree 
Field Galaxy Redshift Survey (2dFGRS) if they re-analyze 
their data using the shallower isophotal limits for galaxy 



Figure 2. Our ansatz for the mass dependence of the disc baryon 
fraction fd . This models a factor of two decline in the value of fd 
from high mass to low mass halos. which fit s the median val ue of 
fd at high masses ( Vlo, Ma o & White 1998) and low masses (van 



den Bosch & S waters 2001) 



magnitude employed by these surveys; see Blanton et al 
(200lQor discussion. 



The theoretical velocity function is given by: 
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much larger than observed, which is the discrepancy we seek 
to resolve. 

This cannot be resolved by making fd a function of 
M200, unless the /<j(M2oo) dependence is unrealistically 
drastic. To illustrate this we assume a different functional 
form for fd(M2oo), as shown in fig. ^. Our ansatz for the 
fd(M2oo) relation models a fairly mild factor of two decline 
in the value of fd from high mass to low mass halos, which 
fits t he median value of fd at high masses ( Mo, Mao fc White 



199c) and low masses (van den Bosch & Swaters 2001). If 



where V c is the circular velocity in the flat part of the ro- 
tation curve and /bright is the fraction of halos which form 



stars. For the mass function dN/ dM, we use the Sheth & 
Tormer (1999) mass function, which gives a good fit to the 
mass function seen in the Hubble Volume simulations (see 
Jenkins et al. (2001)). For evaluating the derivative 



log(V c ) ' 

we use the fitting formula of Mo, Mao fc White (1998] ) , which 
gives the circular velocity V c as a function of the virial mass 
of the dark halo (M200), the concentration parameter (c). 



the spi n of the dark halo (A) and fd- As in Navarro, Frcnk 
& While (1997), for a given M200 we evaluate the concentra- 



tion parameter from the typical collapse redshift. Below we 
discuss our choice of fd , A further and evaluate /bright using 
Toomre instability arguments. 

Let us start by assuming /bright = 1, a constant disk 
mass f raction fd — 0.06 (i.e ., about half the nucleosynthesis 
value (Fukugita et al. 1996)), and the disk spin parameter 



to be the mean value fou nd in N-body simulations A mcan = 
0.042 flBullock et al. 200l| ) . In figure | we plot the theoretical 
velocity function (upper solid line labeled by fd = 0.06). At 
low velocities, the predicted number density of galaxies is 



this was the only process at play it would produce only a 
mild change in the theoretical velocity function, as shown 
by the thin solid line labeled by NDG (no dark galaxies) in 
figure 

To include the effects of the Toomre instability cri- 
terion we proceed as follows. We define a critical spin 
parameter A cr it such that the disc is marginally stable: 
Q[fd(M 2 oo), Acrit, M200, c(M 2 oo)] = 1.5, where Q is evalu- 
ated at one disk scale- length, and we assume a gas sound 
speed c s — 6kms _1 (the gas will in general be photoionized 
by the meta-galactic ionizing background to T ~ 10 4 K). The 
fraction of bright galaxies is given by /bright = f^""* p(X)d\, 
where p(A) is the prob ability distribution of A derived from 
numerical simulations (Bullock et al. 2001). Note that since 
halos with A > A cr it are assumed not to host visible galax- 
ies, the probability distribution of A in observed galax- 
ies is skewed. We evaluate the circular velocity assuming 
a median value A mc dian, defined by the implicit equation 
r median p(A)dA//bright = 0.5. The velocity function we ob- 
tain from computing /bright in this manner is shown in figure 
|l]by the thick solid line labeled by DG (dark galaxies). Thus, 
once we include the suppression of high-spin galaxies by ex- 
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plicitly computing /bright , the theoretical predictions match 
the observations well. 



Th e suppression of accretion after reionization ( Bullock 



2001 



somerville 2001 



Benson et al 



et al. 2001| phiu et al 

2001ah, is an attractive mechanism for regulating fd since 
(unlike for instance, supernova explosions) it affects dark 
galaxies as well. We model this effect by using the fitting 
formula (iGnedin 20001; ISomerville 200 lh: 



ff ot °(M 200 ,z f ) = 



Mm [1 + 0.26 



M C | 
M20O ^ 



(3) 



where we approximate Mc as the mass c orresponding to a 
halo with virial velocity V c = 50kms _1 (somerville 2001) 
at the redshift of formation Zf. Fo r a given halo mass M200, 
extended Press-Schechter theory (Lacey & Cole 1994) gives 
the probability distribution of collapse redshifts p(Zf); we 
can then use equation (^) to obtain the probability distri- 
bution p(fd)- Note that fd = fifc/fim for z > -Zrcion, where 
Zrcion is the redshift of reionization; we assume Zrcion ~ 7 
(the results depend only very weakly on z Icion ). 

In the photoionization model, the fraction of dark galax- 
ies can be computed via equation ^, with the probability 
distribution p(fd) given above. Significant suppression of the 
number density of low circular velocity galaxies only takes 
place off the scales plotted in fig Photoionization only 
successfully suppresses gas accretion in very small halos; it 
succeeds in alleviating the dwarf satellite problem, but leaves 
unaffected halos of V c > 75 Km s~ . By itself it therefore 
cannot produce the low values of fd required to reconcile the 
theoretical and observed velocity function. 

We do not attempt to model the fd(M2oo) relation in 
more detail (which is the task of hydrodynamic galaxy for- 
mation modelling) , but merely note that it is consistent with 
observations and theoretical prejudice that gas expulsion 
and suppression of accretion are more important for shal- 
low potential wells. Possible mechanisms include suppres- 
sion of accretion due to entropy injection in the intergalactic 
medium by supernovae and/or AGN jets, or ram pressure 
stripping as a smaller halo falls into a larger halo. We note 
merely that by invoking the Toomre instability criterion, we 
do not require the very strong variation of fd with M200 re- 
quired in most semi-analytic models of galaxy formation. In 
particular, we do not require a very sharp drop in fd at the 
faint end. 



2.2 The fraction of dark galaxies 

In figure [3] we plot the fraction of dark galaxies /dark as a 
function of halo mass. 

The dashed line is / da rk = 1 - /bright , where /bright 
has been computed in the previous section. The dotted line 
shows the resulting /dark if photoionization was the only 
process at play. As already noted this mechanism only sup- 
presses accretion in halos of low circular velocity and cannot 
reconcile the observed and theoretical velocity function. As- 
suming fd as a deterministic function of M200 might be an 
oversim plification: in reality f d is likely to be highly s tochas- 
tic (e.g., van den Bosch (2001); 



Benson et al. (2001b)). Also 



analysis of rotation curves (JVO02) confirms that there is a 
significant scatter in the relation (see §3). We thus compute 
also the expected fraction of dark galaxies by neglecting any 
possible mass dependence and using the distribution of fd 
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Figure 3. Fraction of dark galaxies /dark as a function of halo 
mass. The dashed line is /dark = 1— /bright! where /bright has been 
computed in §2, by considering the Toomre instability criterion 
and a mild mass dependence of the disk baryon fraction fd shown 
in fig ^. The dotted line shows the resulting /dark if suppression 
of accretion due to reionization determines fd- The solid line is 
the dark fraction derived using the distribution of f d empirically 
obtained by fitting the rotation curves. 



(p(fd)) empirically obtained by fitting the rotation curves 
for the NFW profile (see JVO02 and §3): 



fd c 



dfd P (fd)p(X) 



(4) 



where A cr it and fd crit are defined by 
<2(M 2 oo, fd crit, A C rit, c) = 1. The result is shown in 
figure ^ as a solid line. The two approaches (the one 
described in §2.1 or using the empirically obtained fd 
distribution) yield very similar values for /dark, indication 
of the fact that the main mechanism to create dark galaxies 
is the Toomre criterion. 



2.3 Dwarf satellites 

From fig. ^ it is clear that star formation within dark matter 
halos with masses below 10 10 M© will be almost completely 
suppressed (not visible because they will be Toomre stable) . 
Thus this provides a mechanism to hide these low mass ha- 
los. To make a more quantitative estimate of the viability of 
this effect, in figure w we compare the observed cumulative 
mass function of dwarf satellite with predictions from N - 
body simulations. The solid line is from Moore et al. (1999), 



where the velocity has been converted in ma ss. The gray 



Moore et al 



area i s an estimate of the error from fig. 2 of 
(1999| ). The solid and dashed lines are the predicted cumu- 



lative number of satellites once the probability that some 
galaxies will remain dark has be taken into account. In par- 
ticular we have used the predictions of the solid and dashed 
lines of figure ^ respectively. This prediction agrees nicely 
with observations (triangles). We have not shown results for 
masses below 10 9 Mq because the uncertainties in the pro- 
cedure to obtain the corrected line become too important; 
only a very small fraction (< 1%) of dark matter halos will 
be visible thus making this fraction very dependent on the 
exact value of the predicted number of dark galaxies. 
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Figure 4. Cumulative mass function of dwarf satellites of the 
Milky Way. Small triangles are the observations, the solid line 
with the gray area is the predictions from N-body simulations 
with an estimate of the uncertainty. Finally the solid and dashed 
dashed lines show the predictions once the the probability that 
some halos will remain dark (due to their spin parameter being 
above the critical value) has be taken into account. Solid and 
dashed lines have been computed using the /dark fraction of fig. 
4. 



3 ROTATION CURVES 

In this section we present some indication of the existence 
of dark galaxies from analysis of a large sample of rotation 
curves. 

We use the same set of rotation curves as in JVO02 
where we determined the best fit disk parameters within the 
context of both a disk within a NFW profile halo and within 
a pseudo- isothermal profile halo (see JVO02 for details). 

The two sets of observed rotation curves are: 



• Courteau catalog ( Courteau 1997 ) which consists of 
optical (H a ) long-slit rotation curves for 300 Sb-Sc UGC 
galaxies. 

• 64 spiral gala xies (Sa-Sd) observed in H a by 
fc Williams (2000|) 



Palunas 



In total we use 364 optical rotation curves (to avoid 
any proble ms with the smearing of the bea m in radio ob- 
servations (van den Bosch & S waters 2001)) corrected for 
inclination and with error bars. For these galaxies the scale 
length of the exponential disk (Rd) is also measured. The 
two models have three free parameters: M200, r c , and fd- In 
principle there is another disk parameter: the modified spin 
parameter A'. If A is the spin parameter of the dark halo, 
then A' = jd/fdX where jd = Jd/J and Jd,J are the total 
angular momentum of the disk and halo respectively. A' can 
be readily obtained using the measured value of Rd and the 
best fit parameters (see JVO02). If angular momentum is 
conserved, then A = A', i.e. jd/ fd = 1- 

In JVO02 we used the sample of 364 galaxies for which 
Rd is measured, we found the best fitting parameters (M200, 
fd, and c) for the above two models and recovered A'. Within 
both models for the majority of the galaxies, a meaningful 
set of fitting parameters was recovered. For the remaining 
galaxies we obtained values for fd that where either implau- 
sibly low (f d ~ 0.001) or may be too high (f d £ 0.2). 



3.1 Fit to the disk models 

As described in detail in JVO02, for each galaxy rotation 
curve and for both models, the best fitting parameters were 
obtained using a standard % 2 minimization and exploring 
the whole likelihood surface. The sample of galaxies studied 
is very non-uniform: some curves have much smaller error- 
bars than others, some curves do not extend to large enough 
radii to show the flattening of the rotation curve, while oth- 
ers show strong evidence of spiral arms and bars in the ro- 
tation pattern, moreover we have not attempted to model 
any bulge nor bar component. Thus it is is important to 
keep in mind that galaxies might be more complicated than 
our model for galaxy rotation curves. However, for most of 
these galaxies, the high-quality of the rotation curves mea- 
surements allow degeneracies among the disk parameters to 
be lifted (JVO02). Here, we are not interested in a detailed 
modeling of individual galaxy dynamics but rather the gen- 
eral statistical trends of the recovered disk parameters from 
the whole sample. 

Fig. |H| shows the distribution of the best fitting disk 
parameters in the fd-X' plane and M2oo-fd plane, for the 
NFW profile (top) and pseudoisothermal profile (bottom) ; 
diamonds correspond to galaxi es from the |Courte au (1997) 
sample and triangles are from Palunas & Williams (2000| ). 
In JVO02 we showed that low surface brightness galaxies 
(LSBs) follow the same trend as the rest of the sample, so 
any observational bias against LSBs does not seem to affect 
this result: the only difference is that LSBs have system- 
atically higher values of A' for a fixed value of fd- In both 
models there is a correlation between A' and fd, albeit with 
a scatter, and a weak anti-correlation between M200 and fd- 

For f d < 0.02, A' values are between 0.01 and 0.03, 
while for fd > 0.02, A' takes all values between 0.02 and 0.4. 



The orig in of this correlation is controversial: van den 



Bosch (2001) argues that when feedback by star formation 
is included in models of galaxy formation, such a correla- 
tion between A' and fd arises natur ally for low ma ss galax- 
ies, where feedback is most efficient; Burkert (2000) showed 
that it can be an effect of disk growth happening inside-out; 
or disks might preferentially lose baryons with high angular 
momentum. On the other hand, Burkert, van den Bosch & 



Swaters (2002) argue that the correlation may be entirely 
due to intrinsic degeneracies due to the fact that from rota- 
tion curves analysis the total dark halo mass is not well con- 
strained. In JVO02 (figure 14) we find that this is the case 
for rotation curves that are not measured at large enough 
radii, but, for most of the galaxies in the sample (70%), the 
high-quality rotation curves allow model degeneracies to be 
lifted. This indicates that the random error in recovered pa- 
rameters is small. On the other hand, the systematic errors 
(especially due to the uncertainty in the dark matter profile) 
may be larger, possibly causing the correlations. 

We argue here that the most significant feature of the 
left panels of fig. 1 is not the correlation itself, but the fact 
that spirals avoid the upper left portion of the X'—fd plane 
and that this is in agreement with the Toomre instability 
criterion. In passing we also note that disk instability might 
make galaxies with disk parameters in the bottom left corner 
of the plot to drop out of the sample: for these parameter 
values discs might be too concentrated and form a spheroidal 
instead of a disc or galaxies might become bulge dominated. 
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Figure 5. Distribution of the best fitting disk parameters in the /^-A' plane and in the fd — M200 plane for the NFW (top) and 
pseudoisothermal profile (bottom). We argue that the Toomre instability criterion correctly predicts the zone of avoidance (gray shaded 
areas) for combinations of disk parameters that involve low M200 an d low fd- Note that there seems to be another small void region in 
the bottom right corner of the /tj-A' plot. This might arise because for these parameter values discs might be too concentrated and form 
a spheroidal instead of a disc. 



Tli3 right panels show an interesting, although weak, 
correlation between M200 and fd'- low mass spirals (M200 < 
10 n Mo) tend to have fd > 0.1 since no galaxies are found 
in the lower left portion of the plot. There can be several 
different effects that module the dependence of f d on the 
halo mass. Linear theory evolution from inflation predicts 
a weak anti correlation between spin parameter and peak 
height ( Heavens fc Peacock 198S| ), and hence mass of the 
halo. The expected anti correlation is, however, very weak. 
On the other hand, cooling is less efficient in massive halos 
(-M200 ii 1O 12 M0), thus fd i s expected to decrease strongly 



(2001b)) at the low mass end. The correlation we find from 



with increasing M200 (e.g., van den Bosch (2001)) at the 



high mass end. However, in small mass halos, reheating of 
cooled gas through feedback and photoionization processes 
are much more efficient than in massive syste ms, thus im- 
plying that fd should increase with M200 (e.g., Benson et al 



rotation curves analysis is in rough agreement with the high- 
mass end trend (see fig. 5 of van den Bosch (2001)). For his 
no feedback model, the baryon fraction decreases for increas- 
ing masses. For the feedback model this is only true for M200 
bigger than ~ 5 x 10 Mq. For lower masses this trend re- 
verses and we do not find evidence for it from our analysis of 
rotation curves. This may be due to the fact that for these 
low masses, only those galaxies with highest fd are visible. 

Here we show that Toomre instability criterion predic- 
tions are in agreement with the zone of avoidance for com- 
binations of disk parameters that involve low M200 and low 
fd or low A. 
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3.2 Stability of discs 



The void areas found in fig. ^ can be explained by resorting 



For example discs in which the self gravity is dominant are 
likely to be unstable to the formation of a bar. He re we use 



disc 



to the 'jToomre instability criterion: for Q > 1 the gas 
is stable and thus slai fuimation is strongly suppiessed. Ex- 

pressions for Q and S for the N FW profile have been pre- 
sented in Jimenez et al. (1997), here we derive the corre- 



the results of Efstathiou, Lake fc Negroponte (1982) and 



Mo 



sponding expressions for the pseudoisothermal profile. 

Recall that in the pseudoisothermal profile the density 
p as a function of the radius r is p(r) = po[l + ( r / rc ) 2 ] -1 > 
where po denotes the central finite density. The disk surface 
density as a function of the radius is E(r) = Eo exp(— r/R d ), 
where E = f d M 20 o/{2TrR 2 d ). Appendix I of JVO02 gives a 
relation between po and the concentration parameter c and 
a fitting formula for R d as a function of the disk parameters: 
M 2 oo, po, r c and A'. 

The shaded areas in figure ^ represent the void regions 
predicted by the Toomre criterion assuming c s = 6 km s -1 
for all values of masses and radii. The dashed line shows how 
the avoidance region would change if c s = 2 km s -1 . Such 



a low value is based on the findings of Corbelli & Salpeter 
(1995^7|who argue that when the quasar UV ionizing back- 
ground is the only source of ionization of the HI gas, as 
should be the case for dark galaxies with no supernova, then 
the sound speed of the gas cannot be higher than 2 km s _1 . 
Measured values of the so und speed in the M ilky Way and in 
nearby galaxies (see e.g., |Kennicutt (198S )) range between 
3 and 10 km s _1 . Note that using a low value for c s is a con- 
servative assumption, since for small values of the velocity 
dispersion the formation of dark galaxy is more difficult. 

More specifically, for a grid in M200, A', f d , we compute 
the surface density of the exponential disc E and whether 
Q > 1. The shaded area in the left panels of fig. |^ are ob- 
tained assuming a fixed mass of 5 x 10 11 Mq, approximately 
in the middle of the mass range considered. The shaded area 
in the right panels is obtained as follows: we assume a log- 
normal distrib ution of spin parame ters with mean 0.042 and 
dispersion 0.5 (Bullock et al. 2001), thus 68% of galaxies are 
expected to have A' > Ag g = 0.033. We then assume that 
a combination of f d and M200 parameters correspond to a 
dark galaxy when for all A' > Ag 8 the disc is Toomre stable. 
The Toomre criterion seems to work better for the NFW 
profile than for the pseudoisothermal. 

However, for the pseudoisothermal profile model there is 
no theoretical prediction for the expected A' distribution and 
the dependence of the concentration parameter c = ifeoo/fc 
on the halo mass. To produce the gray areas in fig. [B] we 
assume the same log-normal distribution for A as predicted 
from CDM simulations. This assumption should be valid 
since JVO02 find that the empirically recovered A' distri- 
bution is well approximated by the CDM prediction. For 
setting the concentration parameter we consider that, when 
a galaxy rotation curve is fitted by the pseudoisothermal 
profile, r c obtained is about 1/10 of that obtained when the 
curve is fitted by the NFW profile (fig. 5 of JVO02). Thus, 
for a given mass, we compute the concentration parame- 
ter as expected in the NFW profile, and convert it to the 
pseudoisothermal one. The scatter around this relation is, 
however, quite large, thus introducing some uncertainty in 
the grey areas so obtained. 

The Toomre criterion is a local instability criterion; 
global instabilities can in principle trigger star formation. 



Mao & White (1998) to conclude that the onset of these in- 
stabilities does not affect the grey areas. 

We conclude that the Toomre criteria at least in the 
context of a NFW profile seems to predict correctly the re- 
gion of the parameter space that galaxies tend to avoid (i.e., 
these galaxies should be dark). 



4 DISCUSSION AND CONCLUSIONS 

In this paper we propose that gas in a large fraction of low 
mass dark matter halos may form Toomre stable disks. Such 
halos would be stable to star formation and therefore remain 
dark. This may potentially explain the discrepancy between 
the predicted and observed number of dwarf satellites in the 
Local Group, as well as the deviation between the predicted 
and the observed faint end slope of the luminosity function. 
We show that model fits to rotation curves are consistent 
with this hypothesis: none of the observed galaxies lie in the 
region of parameter space forbidden by the Toomre stability 
criterion. Such Toomre stable disks may be the origin o f 
damped Lya systems seen at high re dshift (Kauffmann 1996 ; 
Jimenez et al. 19981; |Mo et al. 1998J). 



Most semi-analytic models of galaxy formation achieve 
a reconciliation between the observed and predicted abun- 
dance of low luminosity galaxies by drastically decreasing f d 
for faint galaxies. At present, there is no evidence from rota- 
tion curve modelling that low circular velocity disks are dark 
matter dominated (which would be the case if f d were very 
small). Our model only requires a very gentle decline in f d 
toward low masses. It magnifies the effect of previously pro- 
posed mechanisms (e.g. supernovae feedback, suppression of 
accretion), since Q oc f^ 1 and low f d disks are more likely to 
be Toomre stable; such mechanisms can therefore be 'tuned 
down' to lower levels. This may be relevant to the problem of 
trying to simultaneously fit the luminosity function and the 
Tully-Fisher relation. In our model, the mass-to-light ratio is 
much more stable as a function of halo mass (except in dark 
halos, where it is infinite). Indeed, by construction we fit 
both the Tully-Fisher relation and the luminosity function, 
since both were used to o btain the velocity function. 

Jimenez et al. proposed that disks in high spin 

halos may not be low surface brightness galaxies but instead 
be completely dark. This was at variance with the proposal 
that a large fraction of the 'missing galaxies' lie in low sur- 
face bright ness galaxies which formed preferentia lly in high 
spin halos ( Dalcant on, Spcrgcl fc S ummers 1997 ). Here, we 
have extended the !Tnien'eTet'*aT**(1997| ) model to a LCDM 
scenario and show that it is supported from an analysis of 
rotation curves. If it is true that disks in such high spin halos 
may not be low surface brightness but instead be completely 
dark, then there should be a cutoff in the surface brightness 
distribution of galaxies; as surveys probe increasingly low 
surface brightness galaxies the luminosity density will not 
continue to increase but instead will plateau. 

It is important to consider independent indications for 
the existence of this dark galaxy population and devise pos- 
sible observational tests. For example, Augusto fc Wilkinson 
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(2001) and rrcntham, Moller fc Ramirez-Ruiz (2001) have 



investigated the detectability of dark objects. 

These dark galaxies have very low surface densities of 
HI (below 10 20 cm -2 at one scale length of t he disk), stil l 
high enough to be detected by current surveys (Zwaan 2001) 
if the HI was in its neutral phase. Current observations that 
are sensitive enough to detect such amount of neutral H I 
have been carried out recently by Charlton et al. (200C ) 
Zwaan fc Briggs (2000 ); Zwaan (2001). In particular, p^waar 
(200iyound no significant detection of HI, down to a detec- 
tion limit of 7x 10 6 M©, not associated with visible galaxies. 
Therefore, it seems not likely that large masses of neutral 
HI are harbored in dark galaxies. This should not come as a 
surprise, since, given the low surface densities, the hydrogen 
might be ionized by the extragalactic background. At den- 
sities below 10 19 cm -2 , HI is not able to self-shield from the 
external radiation field and most of it will be in its ionized 
phase. Another possibility is that, due to pressure ram strip- 
ping from the hot gas of the hosting large galaxy, most of 
the HI gas could be r emoved from the dark galaxy (|Quilis ' 
Moore Bower 200?] ). 



A more promising route to detect the presence of dark 
galaxie s is through their gravit at ional lcnsing properties. Re - 
cently |Metcalf fc Zhao (200l|); |Dalal fc Kochanek (200i| ); 



Keeton (2001); Bradac et al. (2001) have discussed the need 



of substructure to explain the relative fluxes of multiple 
lensed systems. They show that the flux ratios observed in 
each lens can only be explained by the presence of substruc- 
ture within a large smooth halo. This population would be 
similar to our dark galaxies. 

Dark galaxies with dark halo masses above 10 10 M Q 
should be directly visible via strong gravitational lensing. 
They will produce splitting angles <y 0.1 arc second, which 
is just below detection with current technology, but will 
be easily observable with future radio interferometers (e.g., 
VLBI) 

Using the metallicity dis tributions of globular clusters, 
Cote, West & Marzke (2001) conclude that the mass spec- 



trum of proto-galactic fragments for the galaxies in their 
sample has a power law with index ~ —2, indistinguishable 
from that predicted from N-body simulations. They argue 
that the missing satellite population must therefore belong 
to a class of dark galaxies, similar to the ones we consider 
here. 

The biggest caveat in our model is the assumption of 
conservation of angular momentum. In principle angular mo- 
mentum loss is possible due to torquing by substructure in 
the dark matter halo; indeed, it is invariably seen in SPH 
simulations of disk galaxy formation. This results in more 
compact disks which are much more susceptible to gravi- 
tational i nstability. However, both fits to the scale lengths 
of disks ( [Mo, Mao fc White 1998j) and modeling of the ro- 



tation curves (Uimenez, Verde & Oh 2002 



van den Bosch 



Burken & Swaters 2001) recovers spin parameters consis- 



tent with the log-normal distribution seen in dissipationless 
N-body simulations, indicating little or no loss of angular 
momentum. In addition, we have ignored the effect of merg- 
ers; as halos merge their disks are likely to be disrupted 
and transformed into spheroidal systems, which again are 
more likely to become self-gravitating. The merger history 
of halos is well handled in semi-analytic models of galaxy 
formation which employ Press-Schechter based Monte-Carlo 



merger trees. It would be very interesting to incorporate 
the additional physics of Toomre stability into pre-existing 
semi-analytic (or hydrodynamic) galaxy formation models 
to compare its importance against various other proposed 
schemes for suppressing star formation in low mass halos. 
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